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Recently, the STAR collaboration at RHIC has presented experimental evidence for the correlation 
between the elliptic flow difference of charged pions and charge asymmetry as a possible signal of 
the chiral magnetic wave. We demonstrate that the STAR results can be understood within the 
standard viscous hydrodynamics. 


Recently, the STAR collaboration at RHIC has measured the difference in the elliptic flow parameter V2 between 
7r+ and 7r“ on an event-by-event basis, and found a linear dependence on the charge asymmetry of the collision system 
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where N± is the multiplicity of positively (negatively) charged particles. This quantity has attracted much interest 
over the past several years as a possible signal of the so-called chiral magnetic wave (CMW) Q. It has been argued that 
a strong magnetic field created in off-center heavy-ion collisions together with the chiral anomaly induce a quadrupole 
deformation of electric charges, resulting in the elliptic flow difference ©■ The main predictions of the CMW are 
that AuJ depends linearly on Ach, and the slope parameter r is positive. Moreover, r has a characteristic peak as a 
function of the centrality. The STAR data are qualitatively similar to these expectations. 

Of course, one has to thoroughly examine all the other ‘non-exotic’ mechanisms that can contribute to the difference 
© and subtract their contributions before finally claiming the discovery of the CMW in heavy-ion collisions. Some 
attempts in this direction have been made in [1, Q i but these alternative scenarios are already disfavored by the data 
[H (see, however, [|). In this paper, we raise the possibility that the STAR data can actually be understood within 
the standard viscous hydrodynamics. 

In a previous paper Q , we have analytically computed the difference AuJ for the anisotropic Gubser flow Q coupled 
with conserved currents and found that it is proportional to the shear viscosity 77 and the isospin chemical potential 
/r/. The result is, neglecting numerically small corrections. 
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where K oc rj/s is the Knudsen number (s is the entropy density) and T is the freezeout temperature. (The sign 
convention of Av2 here is different from [^.) In heavy-ion collisions, the mean value of pLj is slightly negative because 
the colliding nuclei are neutron-rich. In place of fj-iy we may alternatively use the charge chemical potential t^Q- 
Numerically, yti are fiq are close (cf. Ref. [8|), and we expect that the following discussion will be qualitatively similar 
in the two cases. 

In order to establish the connection between ©, and we classify events according to the value of Ach and assign 

effective freezeout parameters in each bin of Ac^El This can be done by using the statistical model of hadrons in 
which the multiplicity Ni of hadron species i is computed as Q 


° g (^) exp , (3) 

where V is the volume, rrii is the mass and gi is the degeneracy factor^ The sign =pl corresponds to fermions/bosons. 
B,I,S are the baryon, isospin and strangeness quantum numbers, respectively, and Rb,i,s are the corresponding 
chemical potentials. The asymmetry Ach can be evaluated by summing over all charged hadrons N± = whose 


^ To assign freezeout parameters for certain subevents is a common practice in heavy-ion collisions. For instance, T and p. are often 
plotted as a function of centrality at fixed energy. Experimentally, centrality is determined by multiplicity, and the idea here is that 
events with the same multiplicity are regarded as a statistical ensemble and chemical potentials can be independently assigned for this 
ensemble. Here we do the same, using Ach bins instead of multiplicity bins. In the STAR measurements, the number of events in each 
bin of Ach fixed centrality is typically O(IO^). 

^ The temperature T here is the chemical freezeout temperature which in general differs from the kinetic freezeout temperature in lO. 
Here we assume that the two temperatures are equal or close to each other, following the early freezeout model @ in which the formula 
([ 2 I) was derived. Note that in the same model the ratio /ij/T in |[2j| is approximately constant during the time evolution. 
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FIG. 1. Charge asymmetry Ach as a function of the isospin chemical potential /r/. Blue solid line: T = 159 MeV with the 
parameterization ((6]). Blue dashed line: T = 159 MeV with ([5]). Red line: T = 90 MeV with ([6)). 


masses are below 2 GeV0 This determines a function Achif^B, l^i, /ts, T) which is to a very good approximation lineai@ 
in /r’s, 

Ach^ c{T)^ib + c'{T)^ii + c''(T)fis , (4) 

with all the coefficients (‘susceptibilities’) positive c' > c" > c > 00 Naively, there seems to be a sign mismatch if one 
uses O to rewrite © in the form o because r > 0 and c' > 0. However, one should take into account the fact that 
in heavy-ion collisions /is,/!/ and fis are not independent of each other. The ratios fJ-i/fJ-B and ^-s /^-b are more or 
less universal as they are determined from the quantum numbers of the colliding nuclei, namely, isospin asymmetry 
and strangeness-free conditions Bm. These conditions can be and have been used at different energies, different 
centrality bins and rapidity windows (i.e., different types of subevents). We thus assume that these ratios are fixed. 
More precisely, in practice, we consider the following two parameterizations which we extracted from the result of the 
statistical model fits of experimental data from the SPS to RHIC, 7.6 GeV < ^/snn < 200 GeV [l^ 


0.0308/iB-b2.77- 10 “Vb, 

US = 0.249^3 - 1.09 • 10"Vb , 

(5) 

fii = -0.293- 0.0264^3, 

US = 1.032-f 0.232^3, 

(6) 


where fAs are in units of MeV. The first choice ([5]) which includes cubic terms is motivated by a recent lattice study 
0. The second choice allows for possible intercepts in the limit —>■ 0. Actually, there is some arbitrariness in 

dehning when hb < 0. Yet, keeping the same ratios ni/fiB and hs/ab is natural from the viewpoint of the charge 
conjugation symmetry. 

With the above parameterizations, A^h becomes a function of and T. This is shown in Fig. [T] at T = 159 
MeV relevant to ^/snn = 200 GeV at RHIC [ 13 . We observe that, somewhat counterintuitively, Ach is a decreasing 
function of /i/. Actually, the slope is sensitive to the temperature. As also shown in Fig. [1] if we artificially lower the 
temperature to, say, T = 90 MeV, Ach becomes an increasing function of fij. This indicates that, at high temperature, 
the charge asymmetry is not dominated by pions. 

We thus Hnd that o and ([2]) are consistent including the sign of r. The magnitude of r, on the other hand, 
depends both on the centrality and the collision energy as observed by the STAR collaboration [I|. We confront 


^ In practice, we include the k = 1, 2, 3 terms in the sum m for pions, k = 1,2 terms for kaons and only the k = 1 term for all the other 
hadrons. We have checked that other higher k terms are negligible. 

^ At low energy where ^ 200 MeV, we start to see weak deviations from the linear fit 0 which however do not obstruct the 
determination of r. 

^ At T = 159 MeV, and when the linear approximation is valid for Ach'> we find c = 0.383, c' = 5.45 and c" = 1.45 in units of GeV“^. 
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FIG. 2. Left figure: AuJ versus A^h with T — 159 MeV at ^/snn = 200 GeV, 30-40% centrality. The solid and dashed lines 
correspond to © and ©, respectively. Right figure: The slope r as a function of centrality at -^snn = 200 GeV. The data 
point r = 0.032 at 35% centrality is the input for the theory curve. The STAR data labelled with V2{2} and ii2{4} are obtained 
by using the two-particle and four-particle cumulant methods, respectively. Only statistical error bars are shown in the above 
plots. 


these data with our theoretical results in For this purpose, we first fix the normalization of in ([2]) by 

fitting the slope r = 0.032 measured in the 30-40% centrality range at ^/snn = 200 GeV. We find = 0.28 

and = 0.32 in the two cases ([5]) and (|6]), respectively, see Fig. [2](left). Note that in this fit the ‘intercept’ 

AvJ(base) in ( [p. w hich was experimentally found to be nonzero and positive, is automatically reproduced in the case 
®. See Ref. [m fllj for CMW-based derivations of Ar! 2 (base). 

In the following, we shall focus on the case Let us check whether the number = 0.32 is reasonable. 

For 30-40% centrality, V 2 ~ 0.07 according to [Q. This is related to as (see Eq. (61) of [T^ l 

V2 = (l - , (7) 

which gives K m 2.78. On the other hand, K has been independently estimated in Ref. H- After correcting the 
difference in the normalization, we find Khere ~ lOTArRef.lia]! so that 7fRef.[i3] ~ 0.26. This is indeed consistent with 
the range 1/7 < AfRef.[i 3 ] < 0.5 at = 200 GeV considered in M- 

We now discuss the dependence of r on centrality and energy using the the following formula which can be derived 
straightforwardly from the results of @ 


AfJ ^i/T tj e dN 
Ach ^ Ach s dy ' 


( 8 ) 


where S is the overlapping area of the colliding nuclei and e is the eccentricity. We assume that 77 /s is a constant. 
The part that is sensitive to centrality is then 


e dN 


(9) 


The dependence of ® on the impact parameter b at a/sw = 200 GeV can be read off from Table 1 of 0- Converting 
this 6-dependence into the centrality-dependence [l^, we obtain the right figure in Fig. [21 The agreement with the 
data is remarkably good, except in the most central bin (and the most peripheral bin) where r from the experiment 
is negative!! We however note that the systematic uncertainties of the most central and peripheral bins (not included 
in the error bars in the plot) are very large. In fact, our curve is not inconsistent with the data in these regions if 
systematic errors are taken into account. 


We have seen above that r can become negative if the freezeout temperature is anomalously low. However, this is an unlikely explanation 
of the sign change because at RHIC the centrality dependence of the temperature is very weak Ea. 
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On the other hand, the energy dependence of r at fixed centrality is governed by the factor 

fXi/T dN 

r ^-. 

Ach dy 


( 10 ) 


To evaluate this, we use the known parametrization dN/dy ~ 0 and compute the ratio at each 

value of the collision energy using the corresponding freezeout parameters 0 . The result is shown in Fig. [^together 
with the STAR data. The experimental error bars are rather large especially at low energies, but the general trend 
is consistent with the theory expectation. We see only a mild energy dependence since the ratio is roughly a 
constant along the freezeout curve. 

In conclusion, we have demonstrated that the STAR result At;J oc Ach can be understood within the standard 
viscous hydrodynamics without invoking the CMW. Our scenario can be tested in future measurements as follows. 
The triangular flow difference AriJ satisfies a similar relation to m and we find AwJ /AtiJ /v 2 ■ This agrees 

with the result of Q, although the mechanisms are different. The kaon V 2 difference is dominated by the strangeness 
chemical potential ns'- U 2 ~ '^ 2 ^ ^ @1- Since A^h oc y,s with a positive proportionality coefficient, the slope 

r for the kaons should be negative, in contrast to the pion case, and the magnitude is expected to be larger. The 
situation is similar for the protons. Note that in these predictions V 2 and V 3 are integrated over 0 < pt < 00 . In 
order to properly test them, a wider pT coverage than the one adopted in [I| (0.15 GeV < pr < 0.5 GeV) is necessary: 
Higher harmonics Vn are dominated by pr ^ nT 0 and the average pr increases with hadron mass {{pt) ~ 0.6 
GeV for kaons at yJs^N = 200 GeV (3). Preliminary data for these observables can be found in 0 , but given the 
narrow coverage in pT they are not sufficient to draw firm conclusions. 
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FIG. 3. Slope parameter r as a function of the collision energy at fixed centrality. 
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